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Introduction 


The methods which are discussed in the following paper have been developed 
with the nuclear many-body problem, formulated as a classical quantum mechanical 
many-body problem with only two-particle interaction and without any assumption 
of the origin of nuclear forces, in mind. The success of the shell model has created 
the very simple picture of the nucleus according to which a nucleon moving through 
the nuclear matter is acted upon by a force from all the other nucleons, represented 
by an average potential. On the other hand the shortness of the range of inter- 
nucleon forces implies that one nucleon is interacting strongly with a very limited 
number of other nucleons at a time. 

Considerations of this kind have led people to different approaches to the many- 
body problem. Brueckner and collaborators [1] starting from formal scattering theory 
have developed a formalism which in a direct way takes account of strong two- 
particle interactions of short range in large systems. However, as pointed out by 
Hugenholtz [2] the method of Brueckner is in essence an approximation method for 
the imperfect Fermi gas, and its applicability to a real nucleus of finite size and 
high density is questionable. The many-body problem with hard sphere interaction 
has been treated by several authors [3] but their main interest has not been in the 
nuclear many-body problem. - 

From the arguments given above it seems plausible that it is advantageous to 
use two-particle orbitals in the description of the entire many-particle system. We 
have tried to develop a consistent mathematical procedure for such a description. 
Our starting point is a variational method for the solution of the n-particle Schro- 
dinger equation. Two approaches characterized by different n-particle trial wave 
functions are discussed and the general equations for the determination of (in a 
certain sense) “best two-particle orbitals’ are derived. Both methods may be re- 
garded as extensions of the Hartree-Fock method. The general equations derived, 
though in principle simple, have a somewhat complicated appearance when written 
out explicitly, and other approximate methods are therefore examined. 

For this purpose it is of interest to examine the problem of two interacting par- 
ticles in a common potential. It is very convenient to choose this common potential 
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as an oscillator, because then the Hamiltonian will be separable if center of mass 
and relative coordinates are introduced. Also more physical arguments favour the 
choice of the oscillator potential (See appendices 1 and 2). 


Discussion of trial wave functions 


We shall assume that the Hamiltonian H of the n-particle system considered 
contains only two-particle interaction derived from a potential v, which is strong 
and has a short range: 


H=Sh,+1/2> dvi. (1) 
1 pe i § 


Here h is the sum of kinetic and potential energy for one particle. Let z be any 
number. It is convenient to write H—z1 as a sum of two-particle operators 


H—z1=) > 9,;(2), (2) 


i<j 
where thus 


0,, (z) = (hi +h,)/(n — 1) —221/(n—1) n+ vj;. (3) 
Let x be a complete set of commuting observables for one particle and denote by 
Lei. is ay [HS (4) 


the corresponding representative of the normalized eigenvector | Z) of (1) belong- 
ing to the eigenvalue H of H. 

The usual procedure is to approximate (4) by products of one-particle functions 
(symmetrized or antisymmetrized) 


ke, «+52 | BD ~T Tf (2). (5) 


The effect of the interaction is considered as a perturbation on the one-particle 
orbitals. The only dependence between two particles which is introduced is that 
resulting from the symmetry requirements. In problems involving strong interac- 
tion this procedure fails because the dependence between two particles will be strong, 
though possibly sometimes of limited range. In general perturbation theory will not 
be able to take account of these ‘‘correlation’’ effects. This is very clearly illustrated 
by the problem of two particles interacting strongly in a common potential (Ap- 
pendices 1, 2). 

It is thus necessary to include the effect of the strong interaction in the form of 
the approximate wave function. This should be done in such a way that the ap- 
proximate wave function is exact in the case n=2. It seems natural to divide up 
the n-particle system in disjoint pairs (and one single particle if n is odd) and write 
the approximate wave function as a (symmetrized or antisymmetrized) product 


’ , , 
Cah, «+5 tn | E> ~ py (at, 22) a (3, 24) -.. Pn (wha, hn) (n even), 
2 


, , , , / / (6) 
Bi ess she | Le Oy ay oe) ee Pn-1(Xn—2,Xn-1) M (an) (n odd). 


2 
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The two-particle functions q; (x’,2’’) are assumed to have the correct symmetry. 
As a first approximation they may e.g. be taken as the eigen functions of some 
definite two-particle Hamiltonian 


a(h, +h,)+ v5 (7) 


with numerical coefficients « and f. 
The trial function (6) is by far not the only one containing two-particle depend- 
ence. One of the simplest choices is the following 


<x1,-..,%n|L>~TIf (ai,2;) (l<i<j<n) (8) 
<j 


where the product is over all possible n(n —1)/2 pairs and the two-particle function 
f(x’, x’) is symmetric or antisymmetric. The correct symmetry of the whole product 
is then automatically guaranteed. This is obvious in the symmetric case. In the 
antisymmetric case consider the interchange of particles i and i+ 1. The only change 
occurs in the factor f(x}, 2j.;) the sign of which is reversed. As any permutation 
may be built up from transpositions of neighbouring elements the statement is 
proved. At first sight the trial function (8) seems to be very restricted as it contains 
only one two-particle function. That this restriction is not serious will be shown by 
an example. 

We shall refer to (6) as the disjoint pair function and to (8) as the overlapping 
pair function. 


The variational procedure 
The eigenvalue problem 
H|#>=£|E)y (9) 
is equivalent to the variational problem 
to find the extremum of <#|H|H#)—H(<H|H>—1) over| Hand £. (10) 


Disjoint pair function. In the variational problem we replace |) by the ex- 
pression (6). The variations of the functions ym (x;, x;) are considered as independent 
and subject to no subsidiary condition beside the requirement of correct symmetry. 
(The procedure is thus different e.g. from the usual way of deriving Fock’s equa- 
tions where orthonormality of the one-particle trial wave functions is prescribed.) 
The expression to be varied is now with selfevident notations 


ie Geile os 2 2.0:;(#)S [a> s--| n2vn> (2 even) (11) 


where S is an operator giving correct symmetry to the ket to the right of it. We 
note that 5 5.0,;(#) is a symmetric operator and thus commutes with S. Variation 
i<j 


of the two-particle orbital |!) gives the equation 
of]. - <n en | > > 01,(Z)S| 12>... [n2n>=0 (nm even). (12) 
t<Jj 


One such equation is obtained for every two-particle orbital. These equations are 
all linear in one orbital and bilinear in all the others. The case with n odd is treated 


completely analogously. 
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The equations obtained above correspond to Fock’s equations obtained by an 
analogous one-particle orbit method. The expression to be varied in that case is 


Ci1¢3] +. <n] 220s (B)S|t> --- [n> (13) 
and the equations obtained 
<3] --- <n] 201;(Z)S|1> ..- | > =0. (14) 


By a linear transformation in the n-dimensional space of the one-particle orbits the 
equations (14) are brought on the customary form (orthonormalization) and this 
transformation does not affect the normalized state vector of the whole system 


<410:2(0) {i> 14> —12>1>} = €l >. (15) 


Returning to our equations (12) we note that it is possible to choose here the 
two-particle orbits orthonormalized. It is, however, questionable if a further linear 
transformation of the orbits would in general lead to simpler equations. In prin- 
ciple it is possible to determine the set of solutions to the equations by an iterative 
process analogue to Hartree’s method. However, already in the case of few particles 
this will be very laborious and it is necessary to choose other simpler methods. 

Overlapping pair function. In the variational problem (10) we replace | #»> by 
the expression (8). The expression to be varied is now 


JJJJ ed xh xf ch! )) TTF eis 05)} LSS < at, 05 0 CB) |i’, 0’ >}ATT F(x’ ay')}. (16) 


The two-particle function f (xj, x) is varied without any other subsidiary condition 
than the requirement of correct symmetry. The equation for the determination of 
f (vi, x;) is essentially nonlinear. 


Simplified methods 


The equations obtained from the variational procedure are even for a small num- 
ber of particles so complicated that an exact solution would be extremely laborious. 
As the interaction is strong it is doubtful if a perturbation method or an iterative 
procedure is suitable. For simplification we restrict the admitted two-particle orbi- 
tals in the variational procedure. The first step in this simplification is to choose 
only two-particle wave functions which can be written as products of functions de- 
pending respectively on the sum and the difference of the space-coordinates of two 
Beet (in the following x denotes the three-dimensional position vector of one 
particle): 


: c a4 us oe “f DA A 
fm (228) = Om (7 5) an (7 =r M7 

y2 V2 “a 
In the symmetric case d,, (z) must be an even function, in the antisymmetric case 
dm (z) must be an odd function. The dm (z):3 and 8m (z):s must not necessarily be all 


different. E.g. we can choose all the d,, (z):s equal. 


168 


ARKIV FOR FYSIK. Bd 13 nr 14 


In order to simplify the equation we choose for sm (z) and d,, (z) suitable functions 
depending on a small number of parameters. The functions sp (z) are required to 
have the property that for certain values of the parameters the wave function of 
the whole system describes noninteracting particles in a common potential. The 
deviation of the parameter values from these “shell model’’ values thus describes 
the “correlation” effect of the interaction. To fulfill these requirements we take for 
fm (21, £2) eigen solutions of the problem of two interacting particles in a common 
oscillator potential. The requirement of separability uniquely determines the com- 
mon potential as an oscillator. The constants of the corresponding Hamiltonian are 
chosen as variational parameters. If the interaction is taken to be a square well the 
two-particle wave functions may be written as closed expressions involving Whit- 
taker’s functions (See appendix 1). For a first preliminary discussion we choose both 
dm (z) and 8 (z) as oscillator wave functions but with different parameters. 


‘Discussion of trial wave functions for one-dimensional systems 


In order to investigate the consequences of the discussion given above we con- 
sider the simplest possible example, a system of one-dimensional interacting parti- 
cles, with coordinates 2,,...,%,. Our trial wave functions are chosen of the form 


Pm (B43) = Hg ((2 Vw (#1 +22): 
Eng ((2 Vko (x, —a_)) exp { —k@ (2, + 24)? — ke (x, —4q)?} (18) 


where H,,q (z) and H,,6 (z) are Hermite polynomials. H,, (z) is even or odd de- 
pending on the required symmetry or antisymmetry of Y» (21,2). If ke =ko, 
Pm (#1, Hy) is a wave function for two noninteracting particles in a common oscillator 


potential. 
Disjoint pair trial function. For the ground state of a system of bosons we take 


P (%1, Xp) = exp { — ke @ (a + 2)? — ko (#1 — Xy)?}, 


(19) 
Oe eens ss .» %n)=Nexp{—ko 2 2(a +H) ke 22 (ay — 2;)}. 
For the ground state of a system of fermions (even number n of particles) we take 
Ma tMoe=4m—3(m=1,..., n/2) and me odd. Ifkg—ke =0 the wave func- 
tion is 
E Hy (ay) - «+ Hy (A) 
dives - &n)=Nexp{—(k@ tke) > ai} OE ar ae (20) 
‘i 97 5 (4) ool gas (2a) 


(See appendix 3). As was required this wave function describes n noninteracting 
particles in a common oscillator potential. 

Overlapping pair trial function. For the ground state of a system of bosons the 
wave function is the same as in the previous case (19). For the ground state of a 
system of fermions we take mg =0, mo = 1 ie. 

f(y, ®p) = (#7 — %q) exp { — hq (% + 9)? — ko (a — y)?}. (21) 
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In appendix 3, it is shown that 


Hole) 4: atd awe) 
ie (xj — x) z= | tape towels wate? Cle ep lts SWE R|\6 (22) 
1<i<j<n Htan) te len yen) 


The total wave function is then 


W (a, 22. an)=N I], flena)= 


1<i<j<n 
ak THe (t4) skint Waa) 
i a Hy (%n) «-- Hy-1(Xn) 


If kg =ke this is just the wave function (20) describing noninteracting particles. 

It should be remarked that the discussion given above is independent of the form 
of the polynomials. It is possible to choose instead of Hermite polynomials simple 
powers or any other set of polynomials of the same degree. 

Determination of the ground state of a system of bosons. As indicated above the 
trial wave function is given by (19). In order to make the calculations simple we 
assume that the interaction between the particles is derived from a Gauszian po- 
tential V,exp —f (a;—2;)? and that the common potential is an oscillator. The 
Hamiltonian of the system is then 


h2 a Mm 1 
i Set? 25) Ve —B(2j4-2,)* ¢ 
om rae +2” Dats VGr re Z (24) 
We write the wave function as 
SE eee TYP 
isp eur ypan,) SN ent Oar (25) 


where 


K=4(n—l)(kg@t+ko) and k=4 (kg —ko). 


By elementary integrations the expected value of the Hamiltonian (24) is deter- 
mined as a function of k and K: 


(H>=T+ Vier Vints 


fhe ae 
San dpe 
1 n(K + (n—2)k) 
Voge cmt 
Peed Rae) Eye 
K—k 
K-k+p 


(26) 


1 
Vint = 9 Von(n— 1) 


K and k are to be determined from the minimizing of <H >. It is thus possible to 
represent the expected energy as function of the strength of the interaction (given 
by V,) and the range of the interaction (given by 1/8). 
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If there is no interaction we get k=0 and K=2mo/h. 
As a measure of the deviation of the wave function (25) from the corresponding 
independent particle wave function we take the correlation coefficient 


_ _ ¢%%> 


=, 27 
V<xt> <xd> ss 
which may be expressed as a function of K and k 
—1)k (n— 1) (ke — ko) 
(n— 1) ice (28) 


c" K+(m-2)k (2 n—3) (kg +h) 


From the definition (27) follows that —1<o< +1. In an independent particle de- 
scription 9 =0. |o|=1 means that x, and x, are linearly dependent. 

The extension of the methods discussed in this paragraph to three dimensions is 
straightforward in the symmetric case but the antisymmetric case is more involved 
because it is difficult to find a suitable two-particle wave function from which the 
total wave function may be constructed. The requirement which is hard to satisfy 
is that for certain values of the variational parameters the wave function must de- 
scribe independent particles. We hope to be able to return to this point. 


APPENDIX 1 


Two interacting particles in a common potential 


The system consists of two identical particles with mass m and coordinates q, and 
q, in a common oscillator potential 


5 a o:o-g)= 5 (lai + og) +02 @) (A 1.1) 


where 
O= 0,224 Wy G9+0222 


The interaction between the particles is determined by a potential V (q,—q,). The 
Hamiltonian H of the system is then 


= (1/2 m) (pi + pa) + (m/2) (q,- + HG, + Gy: ++ q,) + V(G,—q,). (A 1.2) 


We introduce the new variables fg and fo by the transformation 


=(1/V2) (fg -Fo). 
The Hamiltonian H can now be written 
H=-H5+Ho, (A 1.4) 


where 
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Ho = (1/2 m) P@ + (m/2)rg@-w-o-Fe (A 1.5) 
Ho =(1/2m) pS + (m/2)Fg- w+ o-Fo + V (V2¥o). 
The two-particle eigenvalue problem 
H|Z>=2|E> (A 1.6) 
is replaced by the two one-particle eigenvalue problems 


H|Eo>=Fa|l#Zo> 


E=Rgt+Eo. (A 1.7) 
Ho |Zo>=Ho|l#o> . 


The state vector |Eo> must have even or odd parity corresponding to symmetry 
or antisymmetry of the two-particle state |. 


We shall now assume that 
o=ol 


V (V27) )=a( r)+B(r )r2 + y ( r)r yo? (r=|7|). (A 1.8) 


The operators (A 1.5) are now both invariant under rotations. For the radial part 
f(r) of the wave function we obtain in the ©-case the differential equation 


a? 2 72 9 
[gpa E+) + m/AA)y he KS + me 


+ (2m/h?) (Bo ~ x(m)|f)=0, (A 1.9) 


Tn an interval where «, 6 and y are constants this equation may be transformed 
into Whittaker’s differential equation 


a2 
fa —(1/4)+k, wt (1/A~ pe)a*l We) =0. (A 1.10) 
The transformation is 
x=)r 
W=V2arf(r). (A 1.11) 


The parameters in Whittaker’s equation are determined from the relations 


22 = (m2 w?/h2) + (2m B (r)/h2) 
4u?—1/4=1(1+1)+(2my (r)/h2) 
22k=m (Ee — «(r))/h. (A 1.12) 
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APPENDIX 2 


Two interacting one-dimensional particles in a common potential. 
Variational treatment 


The Hamiltonian for two one-dimensional particles, coordinates gq, and qy, with 
Gauszian interaction, in a common oscillator potential, is 


" h2 e2 2 Mm w2 3 -; (45%) : 
H= raladtea)* 3 (fit 2) — Voe (A 2.1) 
Introducing " 
ay = (1/V2) (dy +42) (A 2.2) 
v= (1/V2) (11 — 4); 
we get 
h* [d%e! ie? LEE TP aa) oe 
is eee eer a) oe 
7) 
=H, (05 )+H, (+5) (A 2.3) 
The orthonormalized eigenfunctions ¢, (x) of the operator 
Gt ek 
Q= ee ra (A 2.4) 


are 


=(1/22)4(1/n!)2 ( = ad ae, (A 2.5) 
The y satisfy the equations 
Q dn (x) = (mn +1/2) bn (2). (A 2.6) 
If 
A? = (h/2mo) (A 2.7) 


the n:th eigenfunction of a harmonic oscillator with mass m and frequency w is 
1 
(1/A)?¢, (a/A). (A 2.8) 


: F 7) ; 
We approximate the eigenfunctions of the operator H, («, ss) by oscillator wave 


functions obtained from the minimizing of the expression 


F, (0?) = [$s (x/o) H, (« =) bn (x/a)dx/a (A 2.9) 
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over o2. According to the definition (A 2.5) 


sel’) ~aslels) +172) (6) 


We introduce the notations 


an = | |bn(x) Patdx (A 2.10) 
bn («) = ice (x) Peat" de, (A 2.11) 


and obtain for F’, (o?) 


h? (dn 
2mo?\4 


F,, (62) = (n+ 1/2) +™ 2 a, — Vsb0((2)')}- (A 2.12) 


For the determination of a, and 6b, («) we use the identity 
oo 1 os 
S| bn (x) [20 = (1/20)2(1/(1—#))Be 27 18? (A 2.13) 
0 


1 2 
If both sides are multiplied by x? and by e 2*” respectively and the equality then 
integrated over x we get 


Se" {ld (a) 2a2da = Say fS=a40 +t/(1 — t)?) (A 2.14) 
and 


eS + 00 Bee As oo at 
De I ld (x) 2e 2°" da = > bn (ax) "= (La) B® —2t-+1 +a) 2 (A2.15) 


From (A 2.14 and 15) we obtain a, and b,(«) by derivation with respect to tn 
times and then putting t=0. The first few coefficients are 


(A 2.16) 


An=2n+1 


Finally we obtain for F,, (c?) 


ror es2) us (Q)} 
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Vo 
Qs 
Q CA 
0.5 
Es 
0 


Vo 
EY 


Fig. 1. Dependence of the energy # of two interacting particles in a common oscillator potential 
on the strength of the interaction Vy. Range of the interaction @= Mf h/20mo. 

EY ground state energy for fermions (variational method); 

E4 ground state energy for fermions (first order perturbation theory); 

ES, ground state energy for bosons (variational method); 

ES ground state energy for bosons (first order perturbation theory). 


The statistical dependence of the coordinates of the particles is represented by the correlation 
coefficient @ as a function of V, (see (27) and (28)). 


where 


o2 
— 


and Vy»=havrp. 

It should be remarked that A?=(4/2mq) is the standard deviation of the coor- 
dinate of an oscillator in its groundstate, and a, A? is the standard deviation in the 
n:th excited state. 9 may be called the range of the two-particle interaction. 

It is now possible to calculate the energy of the two-particle system as a function 
of the range (8) and strength (V);=fwv,) of the interaction. The result of such a 
calculation is shown in figure 1, and a comparison with ordinary first order pertur- 


bation theory is made. 
APPENDIX 3 
If Py (z), P, (2), ..., Pn (z) are n+1 polynomials with degree not higher than n, 


and if the coefficient of the highest power of z is one in each polynomial, the follow- 
ing identity holds 
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Pate) wes Pa (2) if the polynomials are of 
different degrees 
Po (2). ++ Pn (22) = Hea (A 3.1) 
ee baile if two polynomials are of 
Ps (Sats) gawd a \en41) 0 the same degree 


The identity is a consequence of elementary determinant rules. : 
Let H, (z) be a polynomial of degree n with the coefficient of z equal to one, 


e.g. a Hermite polynomial. If Q(z, ...,2n+1) is an arbitrary antisymmetric poly- 
nomial of degree mine’ it may always be expressed as 
Hote)” ane tieee 
Q (21, «+ +5 2n41) = const. ; ; = const. I] — 2j). (A 3.2) 
Hy (2n41).-- Hn (Zn +1) 
This follows from the fact that it is always possible to write Q(z, ...,2%n+1) as @ 


sum of a finite number of determinants of the type 
Hin (2,) oc Hin (2,) 


Hi (2n+1) « - - Hin (2n +1) 


with constant coefficients. According to (A 3.1) the only nonvanishing determinant 
of this type is (A 3.2). 
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